Quantum phase diagram of fermion mixtures with population imbalance in 

one-dimensional optical lattices 
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With a recently developed time evolving block decimation (TEBD) algorithm, we numerically 
study the ground state quantum phase diagram of fermi mixtures with attractive inter-species 
interactions loaded in one-dimensional optical lattices. For our study, we adopt a general asymmetric 
Hubbard model (AHM) with species-dependent tunneling rates to incorporate the possibility of 
mass imbalance in the mixtures. We find clear signatures for the existence of a Fulde-Ferrell-Larkin- 
Ovchinnikov (FFLO) phase in this model in the presence of population imbalance. Our simulation 
also reveals that in the presence of mass imbalance, the parameter space for FFLO states shrinks 
or even completely vanishes depending on the strength of the attractive interaction and the degree 
of mass imbalance. 
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The pairing of fermions in multi-component fermi 
mixtures is of fundamental interest not only in con- 
densed matter, but also in atomic, nuclear, and astro- 
physics. When the fermi surfaces of the component 
species are mismatched due to unequal densities, pos- 
sibilities open up for exotic pairing mechanisms such as 
the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) 0, 0] and 
the breached pair (BP) states (EH. 

While these interest- 
ing states of matter remain experimentally elusive, con- 
troversies about their stability have been long-standing. 
Recently, this subject has received renewed intense at- 
tention due to atomic physics experimental studies on 
population-imbalanced ultra-cold fermions [5|- The ad- 
vent of techniques for precisely controlling and detect- 
ing ultra-cold atoms is making it possible to systemati- 
cally study the FFLO and/or the BP states with exotic 
fermionic pairing, which has turned out to be impossible 
to do in conventional solid state superconductors. 

Inspired by the experimental achievements, over the 
past few years a large body of theoretical studies have 
been reported for the exotic pairing states in ultra-cold 
atomic fermion systems with population imbalance 0, 0, 
[IJl • It has been found that both FFLO [1] and BP states 
[7| could exist in 3D systems. Particularly, a consensus 
has developed that the stability region of the 3D FFLO 
states is very narrow, making it very difficult to observe 
experimentally. In contrast, due to the Fermi surface 
nesting it could be much easier to observe the FFLO- 
type states in ID or quasi-lD systems @, Q. In this 
work we present numerically calculated exact quantum 
phase diagrams for the population imbalanced ID fermi 
mixtures with unequal masses for the component species. 

Fermi mixture systems with mass imbalance could be 
characterized by the asymmetric Hubbard model (AHM) 
with species dependent tunneling rates. With equal spin 
populations, a phase diagram for the ID AHM has been 
obtained with the renormalization group technique fiol ] 
and possible spin segregation in this model with repul- 



sive interactions has been investigated with bosonization 
and de nsity matrix renormalization group (DMRG) tech- 
niques [ill ]. Nonetheless, it is not surprising that more 
attention has been given to the systems with both mass 
and population imbalances, in which stabilities of FFLO 
and breached pair (BP) states in 3D systems have been 
extensively studied [121]. Very recently, a quantum Monte 
Carlo (QMC) study on finite-size fermi mixtures in ID 
optical lattices with imbalanced populations and masses 
has been reported[l3j], in which evidence of ID analog 
of FFLO-type states has been found. The QMC simula- 
tion also reveals that when the mass difference is large 
enough, instead of an FFLO-type state, the ground state 
of the system will become an inhomogeneous "collapsed" 
state. 

In this paper, we report a time evolving block deci- 
mation (TEBD) numerical study on fermi mixtures with 
unequal masses and attractive on-site interaction in ID 
optical lattices. Our study is complementary to the QMC 
simulation in Ref. 13], in that we obtain the phase di- 
agrams for such fermi mixture systems. Consistent with 
the QMC results, we find that FFLO states are the only 
possible class of polarized pairing states in such systems. 
As the mass imbalance increases, the parameter space 
for FFLO states shrinks and eventually vanishes com- 
pletely. Since we study the homogeneous system at the 
thermodynamic limit, our phase diagram does not have 
any inhomogeneous collapsed state, in contrast to Ref. 
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For our simulation, we adopt the asymmetric Hubbard 
Hamiltonian to model the Fermi mixtures in ID optical 
lattices with unequal masses: 

- ^2K n 1i + n li) + ^2 S ^( n l* -«Ti)i C 1 ) 

i i 

where c a i (cr =| or |) stands for the annihilation operator 



2 



t { =0.1 0.15 0.3 0.5 0.7 1 




-10 1-10 1-10 1-10 1-10 1-10 1 

FIG. 1: Effective magnetic field vs. chemical potential ground 
state phase diagrams for U = —4. From left to right, the 
six panels correspond to t[ = 0.1, 0.15, 0.3, 0.5, 0.7, and 1, 
respectively. The ground states of the system in the regions 
filled with solid dots are found to be of the FFLO type. In the 
shaded regions around the center of each panel, no population 
imbalance presents in the ground states. From — U/2\ = 
+ to fi c , the dominant/subdominant quasi- long range orders 
changes from CDW/SS to SS/CDW. In the remaining regions, 
at least one species of fermions satisfies (n a i)=G or 1 and the 
system behaves as normal Fermi gas. 

of CT-fermion species at site i, n a i{= c^Cd) represents the 
corresponding particle number operator, t a is the species- 
dependent tunneling rate (through out this paper t-\ is 
set to be the unit of energy), U characterizes the on-site 
inter-species interaction strength, and /i ± S/i gives the 
chemical potential of the t/l-fermion. In analogy with 
real spin systems, Sfj, will be referred as effective magnetic 
field in the following, although our system is spinless with 
the t / I components referring to fermions with different 
masses. The unequal mass effects are incorporated in the 
species-dependent tunneling rates through 

t\ m-t ,. 
— = —■ (2) 

We note that the asymmetric Hubbard model defined 
by Eq. (1) could also be realized in mixtures of same- 
species fermions prepared in two different internal states 
by engineering internal-state-dependent optical lattices 

0- 

To investigate the ground state properties of this 
Hamiltonian, we use an infinite lattice version of the 
TEBD algorithm [l5| , which allows us to study the model 
in the thermodynamic limit. A source of intrinsic numer- 
ical error for TEBD is due to the Trotter-Suzuki expan- 
sion (TSE) used in the decomposition of time evolution 
operator. (In our simulation we choose the 4th order 
symmetric TSE.) Furthermore, the convergence of the 
physical results with TEBD is mainly controlled by a 
cut-off parameter x, which characterizes how well one 



preserves the bipartite entanglement of the system when 
truncating the Hilbert space. In this work, we choose 
X = 60. The convergence has been checked to be good 
enough (within ~ O(10 -4 ) for real space correlation func- 
tions) for our purpose, as compared with \ — 80 and 100 
results. 

To identify relevant (quasi-)phases for our model, we 
calculate the real space spin-spin (5*™), density-density 
(D r ), and pairing (P r ) correlations and their Fourier 
transforms Xk = l/VM^ r=Q X r cos(kr), where M + 1 
is the number of sites involved in the transformation (for 
this work, we choose M = 100) and X stands for S, D 
and P correlations. The real space correlation functions 
arc defined as 

qm _ /mm \ _ /„m\ l m \ 
"r — \ B j a i+r/ \ a i /\ B i+r/j 

D r = {uim+r) - (ni)(n i+r ), (3) 

Pr = (CiT C »i C i+r| C i+rT) 

where the spin operators associated with site i is given 
by s™ = with a and (3 =|, | and a m (m = 

x,y,z) standing for the Pauli matrices. In addition, we 
also calculate (c^c CT i +r ) and its fourier transform, which 
gives the momentum distribution for the cr-fermion. 

In Fig. 1, we present the ground state phase diagram 
(effective magnetic field vs. chemical potential) for Fermi 
mixtures based on the asymmetric Hubbard model at 
U = —4. In the shaded parameter regions around the 
centers of the diagrams, all the lattice sites are partially 
filled by both species of fermions and the effective mag- 
netic field S/i is not strong enough to induce any popula- 
tion imbalance. In these regions, the ground states of the 
system have either charge density wave (CDW) or singlet 
superfluid (SS) as the dominant quasi-long range order, 
depending on the filling factor. As the effective magnetic 
field increases, population imbalance is introduced and 
the system could be brought into the regions filled with 
solid dots as shown in the diagram. These regions are 
of central interest for us, since within them the ground 
state of the system is found to be of the FFLO type. (The 
pairing correlation and momentum distribution charac- 
terizing the FFLO-type states will be presented later in 
Fig. 3 for the sample data points marked in the phase 
diagrams.) In the remaining regions of the diagrams, at 
least one type of fermions has zero or unity occupation 
number at all sites ((n a i)—Q or 1) and the system be- 
haves as normal Fermi gas. The green lines in Fig. 1 
show the phase boundaries of the normal states (N). (We 
note that our so-called normal states could be further 
divided into several different categories according to the 
particle occupation numbers and degree of polarization, 
e.g. vacuum state, fully occupied state, fully polarized 
state, etc. But as they are of limited interest, we do not 
distinguish them in this work.) 

By comparing the diagrams corresponding to different 
tj, we can make the following observations. First, the pa- 
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FIG. 2: Effective magnetic field vs. chemical potential groum 
state phase diagrams for U — —10. From left to right, the si; 
panels correspond to = 0.1, 0.15, 0.3, 0.5, 0.7, and 1, re 
spectively. The ground states of the system in the region 
filled with solid dots are FFLO-type states. In the shaded re- 
gions, the fermions of different species have equal population. 
In the remaining regions, the system behaves as normal Fermi 
gas. 
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FIG. 3: Fourier transformations of the pairing correlation 
functions (Pk) and particle number distributions (N-[k and 
N^) for the partially polarized asymmetric Hubbard model 
at U — —4. Pk, JVffc, and N±k are depicted by the solid, 
dashed, and dash-dotted curves, respectively. Other parame- 
ters {tj,, /x, 5/4: (a) {0.1,-2.2,-1.4}; (b) {0.3,-1.7,1.4}; (c) 
{0.3,-2.3,-1.4}; (d) {0.5,-2.4,-2.1}; (e) {0.5,-2.4,0.7}; 
and (f) {0.7, —1.7, 2.6}. Note that the data points correspond- 
ing to these parameter sets are marked in Fig. 1. 



rameter space supporting the FFLO states shrinks with 
the increase of the imbalance in t a or, equivalently, in 
the fermion masses. This can be understood by consid- 
ering the bandwidth. A larger mass imbalance effectively 
leads to a narrower bandwidth (~ At 2 jU), hence when we 
tune ji or 6/i it is easier for the system to fall into the 
nj = or 1 bands and become normal. A second obser- 
vation is that the symmetry of the =/= 1 diagrams is 
different from that of the £j_ = 1 diagram. When t± ^ 1, 
the diagrams are only symmetric about the center point 
(fi = U/2,5[i = 0), reflecting the fact that the Hamilto- 
nian is invariant under the particle-hole transformation 
combined with the inversion about (/i = U/2,5/1 = 0). 
In the tj_ = 1 case, the diagram is symmetric about the 
two axes /i = U/2 and Sfi — 0, since the Hamiltonian 
now possesses an extra symmetry, namely the invari- 
ance under spin flip combined with the inversion about 
O = U/2,5(i = 0). 

We present in Fig. 2 the ground state phase diagram 
for U = —10. One can see that the key features are 
the same as the U = — 4 case. Namely, FFLO is the only 
partially polarized supcrfluid states in the phase diagram 
and its parameter space shrinks with the increase of mass 
imbalance. Furthermore, the shape and symmetry of the 
diagrams also resemble those of the U = — 4 diagrams 
in Fig. 1. Nonetheless, there are noticeable differences. 



With a stronger on-site attraction, it is harder to break 
the pairs. Hence we have higher critical fields 5(A C in the 
U = — 10 case. Besides, the effects of unequal masses 
become more dramatic when the on-site interaction is 
stronger. For example, in contrast to the U = —4 case, 
when 17 = — 10 and t[ = 0.3 one can barely find the 
FFLO phase and at t± — 0.15 the parameter space for 
FFLO completely disappears for U = —10. 

In Fig. 3, we show the pairing correlation func- 
tions and particle number distributions in the momentum 
space for six sample data points with various {t^Sfj,,^} 
in the FFLO regime of Fig. 1. FFLO pairs are known 
to have non-zero center-of-mass momentum and as its 
signature, the pairing correlation function Pk for FFLO 
states has peaks at non-zero momentum \k-\F ~ %f| with 
k a p = {n„i)-K standing for the fermi momentum for non- 
interacting free fermions. From the solid curves in Fig. 3, 
we can clearly see the peaks of Pk at non-zero momenta 
kp = \k^F — fcj.Fl) indicating the presence of FFLO pair- 
ing. One can also see that the particle number distribu- 
tion functions N^k and N^k drop sharply at k^p and k^F, 
respectively. Another noteworthy point is that the mo- 
mentum distribution function for the "heavier" species 
(}-fermion) in an FFLO state clearly exhibits a dip at 
momentum 2k^p — k^p. 

Considering the potential interests in studying the 
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fermi mixtures with unequal masses and attractive on- 
site inter-species interactions loaded in optical lattices. 
We find that with mass and population imbalance the 
ground state of the system could be either an FFLO state 
or a normal state. With the increase of mass imbalance, 
the parameter space for the FFLO states shrinks. When 
the mass imbalance gets too large, the FFLO states are 
no longer stable. 

We acknowledge helpful comments from A. Feiguin and 
W. V. Liu. This work is supported by ARO-DARPA. 
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FIG. 4: Momentum space pairing and density correlation 
functions for the asymmetric Hubbard model at U = —4 and 
ti/tf = 0.15, which is exactly the mass ratio between 6 Li and 
40 K. 



mixtures of 6 Li and 40 K experimentallyfl6| , we also 
present the momentum space pairing (Pk) and density 
{Dk) correlations for the asymmetric Hubbard model at 
U = — 4 and tjt-^ = 0.15 in Fig. 4. The two upper 
panels of Fig. 4 show the pairing and density correla- 
tions in the case with population imbalance, while the 
lower panels show the case with equal populations. From 
the visibility and height of the peaks in the correlation 
functions, one can tell which kind order is more domi- 
nant. First, we look into the case with equal populations. 
When the filling factor slightly deviates from half-filling, 
charge density wave is the dominant quasi-long range or- 
der. As an example, for Sfi = 1.1 and /i = —1.9, we 
have (rii) = 2k^p — 2fc|F ~ 1.09. One can observe a 
sharp peak for Dk locating at k ~ 0.9tt(= 2tt — 2k a p) 
while Pk shows a much broader and lower peak at k = 
indicating that CDW is more dominant than SS. When 
the filling factor is further away from half-filling, the sin- 
glet superfluid order gradually becomes more dominant. 
(The peak of Dk moves toward k = with its visibility 
and height decreasing, while the peak of Pk remains at 
k = and becomes sharper.) Next, we examine the case 
with population imbalance. From the upper left panel of 
Fig. 4, we can see that Pk shows maxima at non-zero 
momenta, which are verified to be directly given by the 
difference in particle number density, signalling the pres- 
ence of FFLO pairing. The results presented in Fig. 4 
indicate that in order to observe the FFLO state in mix- 
tures of spinless 6 Li and 40 K, one should keep the density 
of 40 K well away from half-filling and consider avoid large 
inter-species interaction. 

In summary, we have presented effective magnetic field 
(Sfj,) vs. chemical potential (//) phase diagrams for ID 
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